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^ ' In a stochastic volatility framework, we find a general pricing equation for the class of payoffs depending on 

the terminal value of a market asset and its final quadratic variation. This allows a pricing tool for European- 
style claims paying off at maturity a joint function of the underlying and its realised volatility/variance. We 
study the solution under different stochastic volatility models, give a formula for the computation of the 
Delta and Gamma of these claims, and introduce some new interesting payoffs that can be priced through 
. this equation. Numerical results are given and compared to those from plain vanilla derivatives. 
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1 Introduction 

> : 

The interest of markets in volatility derivatives has significantly grown since the late nineties. The reason of 
such an endorsement by traders and investors is that volatility-linked products allow trading in equity markets 
fNj ■ without necessarily having to take a position in the underlying, thus avoiding the typical problems associated 
\Q | to delta-hedging. As a result, a large body of theory on the pricing and hedging of these products has emerged. 

Neuberger [30] and Dupire [9] pioneered the classic replication of a variance swap via the log-contract. Since 
then the model-independent approach to pricing volatility derivatives has been widely developed, most notably 
• • ' during the last decade (e.g. Carr and Lee [S])- At the same time methods relying on specific models have been 
. ^ ■ proposed; see for instance Matytsin (55], Elliott [TT] et al, Howison [20] et al, Javaheri [22] et al. 

In the years immediately before the last market crisis, new kinds of volatility-related investment paradigms 
begun to arise. Unlike the pure volatility derivatives already present in the market, these new products are 
based on an equity-linked underlying, and use the future realised volatility (i.e. the averaged standard deviation 
of the log-returns) of this same asset as an adjusting factor. In a broad sense, the purpose of this construction is 
to allow a position in the underlying in a classic sense (options, futures) while at the same reducing the inherent 
volatility risk. The most popular investments of this kind are portfolio management strategies generically 
known as target volatility indices, consisting in a periodical review of the underlying exposure according to the 
realised volatility. The purpose is to reduce the vega effect due to the volatility fluctuations for the duration 
of the portfolio. Contracts written on these indices are then in principle priced "as if" the volatility it exhibits 
was constant. This means not only a safer investment but also a much simplified, Black-Scholes alike, pricing 
framework. Examples of target volatility indices sponsored by investment companies are Russell's Controlled 
Volatility Indices , JPMorgan's Commodity Target Volatility Indices, Morningstar Ultimate Stock-Pickers Target 
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Volatility Indexes, DWS NASDAQ-100 Volatility Target Index and BarCap's Revolver. Also, detailed financial 
studies of target volatility indices performances are beginning to appear (Chew, [5]). 

Given the sophisticated nature of the target volatility indices, a first step towards an analytical valuation 
theory of these instruments entails of course understanding how a "vanilla" joint asset and volatility payoff is 
to be priced. Such a vanilla product is the so-called target volatility option (TVO). A TVO is an European- 
type derivative contract that pays off at maturity a random fraction of a vanilla European Call, times a target 
volatility parameter a representing the investor's expectation of the future realised volatility. The random 
fraction considered is the inverse of the realised volatility of the underlying. As one might expect, it turns 
out that at-the-money, this derivative is approximately priced by a Black Scholes Call option having implied 
volatility a. Introducing a volatility adjustment in the payoff effectively determines the reduction of the price 
paid for. Whether one will benefit from this while cashing in the payoff ultimately depends on the accuracy of 
the volatility prediction a. 

A study of the TVO and related pricing methodologies has begun in Di Graziano and Torricelli, [TU]. We gave 
there three distinct pricing methodologies for the TVO, under a stochastic volatility model where no correlation 
is assumed between the asset and the stochastic volatility driving Brownian components. 

In this paper we develop further some of the ideas of [TU], and we do so in two different directions. In first 
place, we remove the assumption of independence between the asset and the stochastic volatility; secondly, we 
find a general pricing technique valid not only for the TVO payoff, but for any sufficiently regular joint asset 
and volatility payoff. 

The latter contribution, far from being a mere theoretical exercise, may very well yield to sensible real-life 
derivative products. As we shall see, by using a volatility correction, it is possible to modify a European payoff 
/ into a claim /, in such a way that taking a position in / will be less costly, and still produce the same payoff 
as / if some predicted volatility event takes place. If an investor wishes to trade in / and has a strong belief 
about future volatility, she may choose trading in / instead, eventually being better off if her prediction was 
correct. 

This naturally brings in the question of finding a way to price the class of general claims F(x, y) of two 
variables: one representing the asset's terminal value, the other one its realised volatility. In a general stochastic 
volatility framework, we find a partial differential equation giving the time-t price of a contract written jointly 
on an asset as well as on its realised volatility, and solve it by Fourier transform methods. 

More precisely, denoting the quadratic variation of the log-returns of a market asset St by 

h = I v u d u , (1.1) 
Jo 

where v t is the instantaneous variance, we are hereby interested in pricing European-style contingent claims 
maturing at time T and having the form: 

F t = F(S t ,It), (1.2) 

for some function of two variables F. Our fundamental variables will thus be the asset and its quadratic 
variation. In a continuous stochastic volatility model the statistical realised variance and volatility of St are 
approximated respectively by It/T and J It /T, so that the class (|1.2I) is completely equivalent to that of the 
joint asset and realised volatility (variance) claims. 

Although our main result is dependent upon the choice and the calibration of the dynamics for the stochastic 
variance, it is universal in the sense that works with any sufficiently well-behaved such model. The possible 
alternative approach, a parameter-free replication pricing like the one advocated by Carr and Lee in [5], al- 
though certainly not prone to estimation errors, would necessarily rely on the famous formula by Breeden and 
Litzenberger of [2], representing a claim on an asset through a portfolio of European Calls and Puts. In practice, 
the issue with this technique is that the market may not offer a sufficient range or density of traded strikes, 
leading to a truncation or discretisation errors in the formula, especially for long maturities. In such instances 
a model- dependent choice, like the one of this paper, may be preferable. 

Similar versions of our main equation are already present in literature. To our knowledge, Lipton [26J first 
gave the Fokker-Planck equation corresponding to the log-price version of (|2.ip . Fatone et al. [T2] gave a 
solution for such backward equation in the Heston model by means of a Fourier inversion, and then used the 
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arising family of probability densities to obtain the price of pure volatility derivatives. Therefore, in their work 
a double integration is needed for a claim depending on a single state variable. Sepp [21] instead presents (the 
log-transformation of) equation (|2.ip for jump diffusions and solves it with a method similar to ours, but then 
he excludes the price variable from the analysis and reverts to solutions for pure volatility derivatives in the 
Heston model. 

In contrast, we will obtain pricing formulae for claims depending also on the final asset price at expiry for 
any well-behaved stochastic volatility model, while at the same time keeping the integration involved to the 
minimum. 

In the spirit of the systematic study by Heath and Schweizer [18J, great care has been take in emphasizing 
a series of sufficient conditions that make the pricing problem mathematically unambiguous. We do this by 
referring to the classic theory of parabolic equations and SDEs (Friedman [T5], [TB], Feller [T3], Kunita [24j ) 
which we hope to revive in the context of financial PDEs. 

The solution approach proposed is the natural 2-dimensional extension of the method for pricing derivatives 
in a stochastic volatility framework introduced by Lewis in |25| . Our fundamental transform will be taken with 
respect to the quadratic variation J, besides the log-price x. Strikingly, only minor modifications in the final 
formulae of Lewis are needed, which accounts how powerful his method is. Aside, it is noteworthy that the 
several-dimension Fourier transform idea can be in principle applied to pricing of claims depending on other 
kinds of non-traded market factors (e.g. Asian options). 

In section 2 we will define our model and derive the pricing equation, which is solved in section 3 together 
with a derivation for the Greeks. Section 4 shows fundamental transforms for various models and discusses their 
existence. In section 5 some claims of the form F(St, It) are introduced, which are then tested numerically in 
section 6. Technical details are proved in the appendix. 



2 Setting up the equation 

The single asset scenario we are assuming consists of a three factor Ito process X t — (St,Vt,It) describing the 
evolution in time of a risky asset St, its stochastic instantaneous variance Vt and its realised variance I t . A 
constant market risk-free rate r exists, and the asset St continuously pays to its owner a proportional constant 
dividend yield d. Valuations relying on such a stochastic variance model are clearly unique modulo different 
choices of a market price of risk, which we hereafter assume to be fixed. This induces a risk-neutral pricing 
probability measure P, which is the only one relevant in all that will follow. Under such a law the price St will 
therefore exhibit a log-return rate of r — d. 

Let (O, P, J 7 , {J-t}t>o) be a market filtration satisfying the usual conditions, and let W} , W t 2 be two J-f 
adapted Brownian motions having correlation pt. The underlying diffusion X t is assumed to be of the form: 

dS t = (r - d)S t dt + ^/v~ t S t dWl 

dv t = a(t, v t )dt + p(t, v t )dW? (D) 
dlt = Vtdt. 

For our purposes we are interested in the behaviour of this process in a finite time range [0,T]. We assume 
the coefficients ait, x) : R5j_ — > R, /3(t, x) : R 2 ^ — > R + to be locally Lipschitz-continuous in x, uniformly in t; that 
is, for all compact sets K C R+, 3Ck > such that: 

sup \a(t,x) - a{t,y)\+ sup \(3{t, x) - (3(t,y)\ < C K \x - y\, Vx,y £ K. ( LL ) 

0<t<T 0<t<T y ' 

In general, (|LL[) is sufficient to ensure that a unique strong solution to (|D|) exists only up to a random exit 
time of R.3_ . To obtain an everywhere well defined solution we must impose the (in the words of Feller) natural 
boundary conditions: 

Px,t 



sup v s = +oo = F x t[ inf v s = = 0, Vx £ R+ , t £ [0, T). (NB) 

t<s<T J ' \*^ S ^ T / 

We also say that the and +oo boundaries for the variance process must not be attainable. Most of the 
commonly used models for stochastic volatility satisfy (|LL|) but not necessarily (|NB|) for every possible choice of 
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parameters. In practice, failure to meet (|NB[) has to be interpreted as the possibility of volatility explosions or 
volatility vanishing. When this happens multiple solutions for X t are possible. However, while treating specific 
models we will impose suitable conditions under which (|NB|) is satisfied, so that X t will always be unique. 

Now suppose we want to trade a derivative that pays off at the maturity date T a certain function of two 
variables: the underlying terminal asset value and the quadratic variation accumulated over [0,T] . The payoff 
is then represented by the random variable F(St,It), where F{x,y) is an integrable function in the joint 
distribution of St and It for all the possible sets of spot variables. By the usual dynamic hedging argument we 
set up a portfolio that is long the contract, and short certain amounts of the underlying and of another variance 
dependent contract. By choosing the hedge ratios as to cancel the portfolio randomness, we argue that under 
no-arbitrage, for the given market price of risk, the portfolio process must earn the risk-free rate r. The time-t 
value V(St, It,vt,t) of the contract can be thus seen to satisfy the parabolic equation: 

dV , „BV dV dV vS 2 d 2 V (3 2 d 2 V „ dV 

- + {r - d )S- + a- + v- + —^ + ^ +P ^S—-rV = 0. (2.1) 

Pricing a cross asset-quadratic variation derivative F(St, It) therefore amounts to solving the Cauchy free- 
boundary problem (|2.1[) in x [0,T] having terminal condition: 

V(S T ,I T ,v T ,T) = F(S T ,I T ). (2.2) 

This is the generalised version, in the price variable, of equation (13) of Sepp |31] . and the dual of the 
Fokker-Planck equation appearing in Lipton [26\ and Fatone et al. [12J. 

As it happens when dealing with parabolic equations arising from financial modeling, results of exis- 
tence/uniqueness of a solution may not be readily available from the standard theory of parabolic equations. 
Typically, this is because of two reasons: coefficients constraints are not met, or terminal conditions (payoffs) are 
not continuous. However, even if solvability remains an issue one has to live witf0, uniqueness of V essentially 
carries over from the uniqueness of the underlying diffusion X t , which is in turn enforced by assumptions (|LL[) 
and (|NB|) . 

Another interesting element of discussion is whether, and under which conditions, using the Feynman-Kac 
Theorem to link the discounted risk-neutral expectation of the payoff to the pricing equation can be considered 
to be equivalent to the no-arbitrage derivation^. This is a standard requirement in the literature; the following 
Proposition then motivates our assumptions on (D): 

Proposition 1. Under assumptions (LL) and (NB) there exists at most one C 21 solution to problem 112.1]) - 
h2.2\) : if such a solution does exist, for x = (St, v t ,It) it is given by: 

V(x, t) = E x . t [e-^ T ~^F{S T , It)] ■ (2.3) 

Therefore under (|LL|) and (|NB[) pricing a claim of the form ()1.2[) is a well posed problem, provided that (|2.1[) 
is solvable. 

We finally impose a few last growth constraints, this time directly on the solution V: 

V(S) <ifi(l + S" 11 ), V(I) < K 2 (l + I> 1 *) (GC) 

for some K 1 (I),K 2 (S) >Q,h 1 ,h 2 > 0. 

The reason of this assumption is technical in nature and will allow to perform the necessary reductions 
while solving the equation. The classic theory of parabolic equations (Friedman [IS], [TH]) provides sufficient 
conditions on the problem itself under which (|GC|) holds Q16J, Th. 4.3, p. 147). However, it is difficult to give 
a comprehensive set of such assumptions in a financial setting, owing to the lack of the necessary regularity of 
many cases of interest: namely, superlinear growth of the coefficients of (D) or discontinuity of F . Alternatively, 
(|GC[) can be checked case by case, for example by using estimates along the lines of those derived by Bergman 
et al. [3] (Th. 1). In any case, this condition is easily seen to hold for most of the cases accounted in section 5 
(see the appendix). 



1 See for example Andersen and Piterbarg [T] on the non existence of moments in the Heston model. 

2 As Lewis showed in |25] (ch. 9) this is not the case in models in which volatility explosions occurr with positive probability. 
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3 Solution to the PDE 



We shall characterise the solution of the PDE by identifying a fundamental transform for the problem, which 
loosely speaking is nothing but the characteristic function of the fundamental solution of the equation. To do 
this we will apply the Fourier transform to (|2.ip with respect to both variables / and log S. Once a fundamental 
transform has been found, we can invert it on a suitable domain of C 2 and then conclude from Proposition [T] 
the existence of a unique price for F. 

Let V(t, S,v, I) be the solution of (|2.1[) and consider the substitutions: 



r = T-t 

x = logS+(r-d)(T-t) 

e r( - T ~ t '>V(S,I,v,t) if y > 
if y < 0. 



(3.1) 



W(x, y, v, t) = 

Equation (|2.ip can then be seen to be equivalent to the problem 



v fd 2 W dW n dW\ n ^dW dW /3 2 d 2 W dW 

1- 9 -I- 1- rv 1- = 



with initial condition: 



2 \y dx 2 dx dy ) ^ dxdv ^° 9» 2 dv 2 dr ' fi-Q 

. F(e x °,y ) ify o >0 , . 

W(x ,yo,v Q ,0) = i Q . {y ^ Q (3.3) 



For (77, id) € C 2 , let the two-dimensional Fourier transform of W(x, y, v, r) be: 

W(id,n,v,r) = { e ixuJ+im W(x,y,v,T)dxdy. (3.4) 



We denote derivatives by subscripts. Consider the transform W T of the time-to-maturity derivative of W. 
By substituting (|3.2[) in the integral above and integrating by parts we find that 

W T = W T , W x = -iidW, W xx = -id 2 W 

W y = -ir]W, W V = W V , W VV = W VV (3.5) 
W xv = -iidW v . 

provided that e lu,x W{x)\tZ = e lulx W x {x)\t™ = e lvv W(y)\ +00 = holds true for some id, n. These relations are 
clear if we know V to satisfy (|GCj) . which then yields p.5p in a 2-strip Qx = {ai < Im{id) < a 2 , Im(r]) > 0} C 
C 2 . 

Fourier-transforming both sides of (|3.2[) and substituting the above relations we have the fundamental PDE 
for W: 

-(a - iidy^pfi) - - (id 2 - iid + 2ir])W = — — . (3.6) 



2 dv 2 dv y v ^' 2 V (/ 5r 

A fundamental transform H(u),t],v,t) for (|3.2p is a solution to p.6p such that H (id, 77, u, 0) = 1. Assume 
that such a solution exists: is nothing else than the (sign-shifted) characteristic function of the transition 
probability density associated with the process (log St, It), and it is thus a holomorphic functional on a certain 
multi-strip O2 C C 2 . 

Denote the Fourier transform of the payoff in the log-price and quadratic variation by F(tJ, rf) :— W(id, r], v, 0), 
itself a holomorphic function on a third multi-strip O3 C C 2 . Since F((d,rj) does not depend on the variable v we 
see that the product H(id,rj,v,T)F(id,rf) is also a solution to VS. 6')) having initial condition F(id,rf). Therefore 



3 See Lukacs [28], (Th. 7.1.1.). Since we are not confined to real arguments, we need not to consider analytical continuations 
around 0, and may instead develop around any point in whose neighbourhood H is holomorphic. This means that O2 will exists 
somewhere, even if in general it may not contain the real axis. 
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by taking the Fourier inverse of H (u>, T), v, t)F{uj, rj) on a multi-line 

3 

T, = {(uj,r)),uj = s + ik 1 ,r) = t + ik2,s,t€R} cO= f]n h k 1 ,k 2 eR, (3.7) 

i=i 

and finally unwinding the variable change we are led to the solution of (|2.1I) : 

— r(T-t) rikx+co rik 2 +oo 

V(S, I, v, t) = 2 / / S^e-^-^^e-^Hiuj, r],v,T - t)F{uj, r))dudr). (3.8) 

Jiki—oo Jik2 — oo 

Finally, Proposition [T] establishes that (|3.8p is the unique price of F(St,It)- Of course, such an argument 
is meaningful provided that a common domain of holomorphy C C 2 of H and F actually exists. 
We summarise all of the above discussion in the following Proposition: 

Proposition 2. Assume that the solution X t of fD\) is such that the dynamics for Vt satisfy iLL)) . fNB\) . and 
that IIGC\) holds. Further assume that fl ^ and let k\ , k 2 € K be such that: 

S = {{u),rf),u) = s + ikwn = t + ik 2 ,s,t e M.} C Cl. (3.9) 

If a fundamental transform H(u>, rj, v, r) can be found, the price of a claim F written on St and I t is given by 
equation i3.8\) . 

This formula is completely general: in principle, under the given assumptions, it allows pricing under any 
stochastic volatility model. 

Another attractive feature of equation (|3.8p is that it allows us to separate, by means of H and F, the 
pricing information coining from the model from that coming from the payoff. This was one of the original main 
contributions of Lewis's work on stochastic volatility models and is equally valid here. Changing the stochastic 
volatility or the function to be valued only requires changing the corresponding transform to be used in Q3.8[) . 
and not the whole re-computation of the solution. 



3.1 Greeks 

The representation found also allows for a straightforward computation of the Greeks. Calling J(u>,t],v,t) the 
integrand in (|3.8p and differentiating V under integral sign we have that the Delta for the contract F is: 

A= — = — — / / —J(uj,ri,v,T)dudn. (3.10) 



J(oj,r],v,T)du!drj. (3-11) 



9S 47T 2 Jikx-oo Jik 2 -oo & 

Likewise, the Gamma is seen to be given by: 

d 2 V e~ r( - T ~^ r kl+ °° r kz+0 ° iuj - lu 2 



f)2 C 4_2 / / C>2 



The derivative A is one the two coefficients to be used in the hedge ratios yielding equation (|2.ip . The 
sensitivity to the initial inatantaneous varinace dV/ dv can be sometimes expressed in a similar fashion as (|3.10[) 
and (|3.1ip , for example in afRne models. Clearly, the ability to fully hedge will depend also on the possibility 
to identify a fundamental set of securities for the market. However, the situation is no more general than that 
of a standard stochastic volatility set-up, because the diffusion (D) shows no more randomness than a model 
with an asset process and stochastic volatility only. 



4 Model-specific fundamental transforms 

We analyse here in more detail the fundamental transforms of the Heston, the 3/2 and the GARCH models. 
The analytical tractability that characterises these in a standard stochastic volatility scenario carries over when 
realised volatility comes into the picture. Remarkably, the solution of ()3.6p depends on the coefficient of the 
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linear term as a parameter. As the variable r\ appears only in such coefficient, the derivations are formally 
identical to that already present in literature, to which our formulae reduce when r\ — 0. 

Following are the transforms, together with their domain of holomorphy as functions of two complex variables. 
Being a Fourier integral, H is everywhere holomorphic in its domain of definition. Complex square roots are 
always understood to be the positive determination. As per our initial assumption, the parameters for the 
models already incorporate the market price of risk adjustment. For a sketch proof of the derivation of (|4.2|) 
consult the appendix; a complete treatment is to be found in Lewis |25) . 



Heston model 

Parameters are assumed to be constant. In such a case a lot is known about the model (Heston [TH], CIR [7J) 
and a fundamental transform can be easily obtained. 

• Dynamics of the instantaneous variance: 

dv t = K{6-v t )dt + ey/v t dWt, (4.1) 

with K,8,e > 0. Volatility explosions never occur; taking 2k6 > e 2 ensures that the boundary is not 
attainable. Thus under these conditions assumption (NB) is met. 

• Fundamental transform: 

H(u, T], v, t) = exp[C(w, i], t) + vD(u, r), t)] 

C(u, r,, t) = ^ (r(&M - d(u, r,)) - 2 log 1 _ c{u ^ ,V) 

b(uj) +d(ui,i]) 
b{w) - d{u),r]) 
b(cu) = k + ieujp 



d(oj, rj) = \/b(Lu) 2 + ( 2 {lu 2 — iuj + 2irj). 



The expression for C uses the argument by Lord and Kahl [27] to avoid discontinuity issues in the complex 
logarithm. The only singularities occur when 1 - c(uj, Tj)e d ^'^ T = causing divergence in both C and D: 
hence the domain of holomorphy of H is C 2 \ S Kt e,p,r where 

S K< e <P ,T = {(u,v) e C 2 | e - d ^)r = C ( W) v) y (4 , 3) 

3/2 model 

We consider the general form as introduced by Lewis [25 : 
• Dynamics of the instantaneous variance: 

dv t = K{9v t - v 2 t )dt + ev t 3/2 dW t . (4.4) 

Whenever 2k > — e 2 we have that +c>o is unattainable; the boundary is natural for any choice of 
parameters. 
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Fundamental transform: 



/2k9 
\ e 2 v ' 



X(x,t) 



k6t 1 F 1 



2k9 

a, (3, —X ( —^—,k9t 



e* - 1 

a{u3,rj) = c{u),rj) - b{w) 
/3(u},rj) = 1 + 2c(u),rf) 
b(u) = (k + e 2 /2 + iujpe)/e 2 
c(u,rj) = \/&(w) 2 + d{ui,r]) 
d(ui, rj) = 2{u 2 — ioj + 2ir\) j e 



(4.5) 



iFi(a, /3, z) is a confluent hypergeometric series and T is the Euler's Gamma function. Since (3 cannot be 
a negative integer the poles of V are avoided, so that the domain of the transform is the whole C 2 . 



GARCH model 



We only take in account a particular instance of this model, namely when the dynamics are simply those of 
a geometric Brownian motion with drift; we also assume p — 0. The case p ^ can be obtained by a simple 
modification of the derivation in |25) . 



• Dynamics of the instantaneous variance: 

dvt = 9v t dt + evtdWt, 

with e > 0. Clearly condition (|NB[) is always met. 

• Fundamental transform: 



(4.6) 



2 /3+i 



P 



-<8/2] 



{' 3 <°} wri 



3=0 

(3 + iz 



jT(l-/3-j) 



1 

+ 4^ 

P = 29/ e 2 - 1 

d(u,ri) = 2^2(lj 2 - iuj + 2ir/)v/e 



'zsmk{zn)K lz {d{uj : ri))e- {f32+z2) ^ T/ *d2 



(4.7) 



Here 1 is the indicator function and K x the modified Bessel function of second kind. By use of the 
appropriate series representation, K x can be extended to an entire function. So we see that H is a 
holomorphic on C 2 \ {uj 2 — iuj + 2ir/} whenever 29 > e 2 , and it is everywhere analytical on C 2 otherwise. 



5 Some joint asset/volatility derivatives 

We present here a list of European-style derivatives paying off a joint function of a terminal asset value and its 
realised variance or volatility. As explained in the introduction, these can all be considered as volatility-modified 
versions of classic payoffs, where the volatility factor reduces the initial price without effecting the payoff if the 
investor's volatility foresight happens to be correct. The target volatility option is one currently traded product 
of this kind. 



Target volatility option 

A target volatility call option is the option to buy a certain fractional amount of shares if the underlying is worth 
more than the strike price at maturity. Such amount is stochastic and depends upon both the target parameter a 
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set when writing the contract and the volatility realised by the asset in [0, T]. Under the independence hypothesis 
between the Brownian motion driving the underlying and the process for the instantaneous volatility, the value 
of an at-the-money Call TVO is approximately the Black-Scholes price of a call option of constant volatility ef 
(see [TO]). 

• Payoff 

F{S T ,I T )=a x f^{S T -K)+. (5.1) 
V it 

• Payoff transform in log-strike and quadratic variation 

F(u,r)) =a(l+i)J — r , for Im(u) > 1,/mfa) > ■ (5.2) 

y Zl] [UjJ — LO z ) 

As is the case for vanilla options, a target volatility put will have the same payoff transform as a call, but in 
the domain we will instead have Im(w) < 0. 

Double digital call 

A derivative delivering at maturity a unit of cash if both the underlying asset and its realised variance at T are 
above two strike levels K\ and K%. In practice we are adding a further strike threshold to a digital call option, 
so it is intuitively clear that this derivative must be priced less than it. 

• Payoff 

F{St,It) — 1{s t >k 1 ,i t /t>k 2 }- (5-3) 

• Payoff transform in log-strike and quadratic variation 

tsiuj e iTK 2 ri 

F{uj, 77) = i for ImUS) > , Imirj) > 0. (5.4) 

U>7] 

Clearly any other Put/Call combination in the two variables can be imagined. 
Volatilty capped call option 

As in the previous example we can cheapen the price of a European Call by adding the further constraint that 
the payoff is not triggered if the terminal realised volatility is not within an acceptable range. 

• Payoff 

F(S T ,I T ) = (S T -K)+I {K ^ VJ ^ £K2} . (5.5) 

• Payoff transform in log-strike and quadratic variation 

I ^VK.T _ ^K 2 T\ for j, x > x Im , s > Q / 5 _g) 

V / [w + iu} 2 )r] 

A more natural version of this product is obtained by requiring that to get a positive payoff the volatility 
should never leave an interval [K\, K2] at each given time t < T. The resulting derivative is a "volatility version" 
of a double barrier option; pricing it therefore amounts to solve a boundary- valued version of (|2.1j) - (|2.2[l . This 
escapes the pricing frame presented; however, such a payoff could be of interest for future research in a context 
of volatility path-dependent claims. 

Volatility struck call option 

This product gives the writer the option to buy an asset at maturity for a notional amount N, times the realised 
volatility of the underlying. The more the stock is subject to shocks, the less likely is the option to be triggered; 
hence an investor could enter this contract if he is expecting low volatility levels. Just like for a TVO, predictions 
about the future realised volatility a are reflected in setting the notional N. 
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Payoff 

FtSrr. Tt) = I Sit — Nil 

T 

Payoff transform in log-strike and quadratic variation 



F(S T ,/ T )= | St-Wi/^I . (5.7) 



T J \ 2 J iui- uj 2 

for 1 < Im{uj) < 3, 7m(r?) > 0. (5.8) 



6 Numerical testing 



Computations for the payoffs introduced in section [3] have been carried out in a Heston model with parameters 
from Di Graziano and Torricelli [10] . in order to compare the values obtained there with the new methodology. 
The underlying process for the variance is given by 



dv t = k(6 - v t )dt + rj-y/vidWt, (6.1) 

with: 

rc = 0.5, 61 = 0.2, r] = 0.3, u to =0.2. (6.2) 

For different state variables and sets of parameters defining the claims, we compare a MATLAB® Monte 
Carlo simulation based on an Euler scheme with sampling from the log-normal distribution (Broadie and Kaya 
@]), against a MATHEMATICA® implementation of (jSHJ). For the TVO, fi gures from the Laplace transform 
pricing method described in |10) are provided. A comparison of the prices of the products introduced with 
their vanilla counterparts is also given; it is striking how much cheaper the new claims are. Nevertheless, under 
favourable volatility scenarios, they produce the same payoffs as their standard versions. 

Table 1: T = 3, t = 0, a = 0.1, S = 100, r = d = p = 0. TVO valuation for different strikes. 



K 


Laplace 


Monte Carlo 


PDE 


Vanilla 




transform 


simulation 


pricing 


Call 


60 


11.3919 


11.3897 


11.3909 


40.0061 


80 


8.7301 


8.7281 


8.7299 


20.7211 


100 


6.7416 


6.7415 


6.7415 


6.9013 


120 


5.2672 


5.2618 


5.2672 


1.4252 



The Black-Scholes price of a Call option of constant volatility a is given for comparison (see |10| . sec. 7). 



Table 2: T = 5, t = 2.5, r = 0.08, d = 0, a = 0.1, S t = 100, K = 85, h = 0.46. TVO prices for various 
correlations. 



p 


Monte Carlo 


PDE pricing 


Vanilla Call 


-0.8 


10.3154 


10.3975 


41.5145 


-0.4 


9.9415 


9.9505 


41.3683 





9.4398 


9.4549 


41.1688 


0.4 


8.9645 


8.9059 


40.8992 


0.8 


8.3136 


8.3025 


40.5433 



Here we compare to a European Call option in the Heston model having same parameters. 
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Table 3: T = 2.5, t = 1, r = 0.1, d = 0.01, Ki = 100, K 2 = 0.24, S t = 120, p = 0.2. Double digital call for 
different realised variance levels. 



It 


Monte Carlo 


PDE pricing 


Vanilla digital Call 


0.2 


0.0951 


0.0943 


0.5358 


0.3 


0.2366 


0.2426 


0.5358 


0.4 


0.4393 


0.4395 


0.5358 


0.5 


0.5335 


0.5330 


0.5358 



Note how the prices converge to that of a digital Call as K2 becomes more likely to be hit. 

Table 4: t = d = 0, T = 2, r = 0.07, St = 110, K = 100, K ± = 0.2, p = -0.3. Volatility capped call option 
prices for different values of K 2. 



K 2 


Monte Carlo 


PDE pricing 


Vanilla Call 


0.35 


7.7812 


7.7743 


37.2632 


0.4 


16.3226 


16.3006 


37.2632 


0.45 


25.1122 


25.0732 


37.2632 


0.5 


31.6069 


31.5497 


37.2632 



The reference vanilla Call has same parameters as the volatility capped call. As the gap between K\ and K2 widens, the same 
behaviour as in Table 3 is shown. 

Table 5: I t = 0.18, t = 1, r = 0.05, d = 0.02, S t = 50, iV = 150, p = -0.5. Volatility struck call option prices 
for different maturities. 



T 


Monte Carlo 


PDE pricing 


2 


4.8291 


4.8815 


3 


8.9873 


8.9383 


4 


11.8885 


11.9086 


5 


14.1661 


14.2002 



7 Conclusions 

In this paper we have explained possible reasons for the introduction in the markets of derivatives written 
jointly on an underlying asset and its accrued volatility. We have discussed the general problem of pricing 
European claims depending at maturity on such an asset and the total quadratic variation it exhibits. A pricing 
PDE has been derived, and a universal model-dependent solution has been found and characterised in terms 
of the model and the payoff. Issues of uniqueness and existence of such a solution have been addressed and a 
detailed mathematical discussion of the domains of holomorphy of the involved functions has been carried out. 
Furthermore, we have provided an easy analytical representation for the Greeks, and given formulae for specific 
models and payoffs. 

Numerical tests support our main result. In addition, figures confirm the intuition that it is possible to 
conceive volatility modifications of liquid instruments, less expensive than the original product, but paying off 
the same amount in market situations an investor may want to exploit. 

8 Appendix 

Proof of Proposition [T] 

The proof is an adaptation of Heath and Schweizer [18] (Th. 1), the existence of a solution being assumed. See 
also Friedman [IB] (ch. 5-6). 
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Let D n be a family of smooth bounded domains invading , and assume V(x,t) is a C 2,1 (]R+ x [0,T)) 
solution to (|2.ip ~ (|2.2p . Let r„ = {inf t\X t $ D n }. By Ito's formula it can be readily seen that if X t is a solution 
to (P) then for all x€D n ,t<T, 

V n (x,t) =E x>t [e- r ( T -*V(X T Ar B ,T , AT n )] (8.1) 

is a solution to the differential problem (|2.ip with boundary condition V^(x,i) = V(x,t), x G <9-D n , and terminal 
condition V (Xt i\ Tn -,T At„); furthermore, it is the only solution there by the weak maximum principle. By (NB) 
and (LL) we have that the probabilistic family yielding (|8.1j) is strongly Markovian, so conditioning (|2.3p at 
time r n and then taking the expectation shows that V(x, t) coincides with V n on D n for all n. Hence, V satisfies 
(|2.ip — (12. 2p on the whole space. Finally, again (LL) and (NB) imply (see for example Feller [13J, Kunita [24]) 
that U{, hence X t , is (weakly) unique and finite almost surely, and this proves the claim. 

Condition (GC) for the payoffs of section 5 

Let V(S,I) = E X) j [e~' '^ T ~^F(S, 7)1 where F is a payoff from section 5. For a double digital call we have 
F{S, I) < 1 so that V(S, I) < 1 and (|GC|) is trivially verified. If F is a volatility capped call then VI: 

V(S, I) < E x , t [ e - r ( T -') (St - K)+} < S, (8.2) 

where the last inequality follows at once from [3] (Th.l) and its generalisation on page 1600. Also, by arguing 
that the no-crossing property (Lemma p. 1577) for St must also hold for the augmented diffusion (D), we see 
that Theorem 1 can be extended to a payoff F(S, I). In particular if F is given by (|5.7p we see that: 

Hence V{S, I) < S, VS, I. Similarly for a TVO we have: 

9V ^<^VTfl, (8-4) 



dS 

which implies V(S, I) < C\S for some C\ > 0, because V is bounded around 1 = 0. 

Proof of equation ()4.2|) 

We make the ansatz: 

H(w, T), v, t) = exp[C(w, 77, t) + vD(u, 77, r)[ (8.5) 

and substitute this in Q3.6P with parameters from the Heston model. One obtains the decoupled ODEs for C(r) 
and D(t): 

a = k6d 

D' = ^—D 1 + D{k + iepcu) - i(w 2 - iu + 2in). 

The Riccati equation for D is solved by switching to an associated linear second order ODE for its logarithmic 
derivative. C is then found by direct integration. 

Coniditon (NB) for the accounted models 

For the GARCH model the result is trivial. The most convincing way of checking (NB) for the other models 
described is using Feller's explosion test (Feller, |14|). In our case, the scale function for the variance dynamics 
in the Heston model is: 

p(a.) = 1 ^ u -^o/^ e ^(u-i) ^ z 2^/^-i e -^{ z -i) dz ^ du (8 7) 
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A necessary and sufficient condition for the process to attain or +00 is that p(0), p(+oo) < +00. As u — > +00 
the integrand is exponentially divergent, whereas convergence in is happens if and only if 2n8 < e 2 . For the 
3/2 model we have instead: 

p{x) = | jT w 2«/e^(l/ u -l) ^ _ 2K/e2 _3 e -^(l/,-l)^ du ( g g) 

so that this time it is p(0) = +00 for any choice of parameters, and p(+oo) < +00 if and only if 2k < —e 2 . 
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